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Abstract 

In a recent letter (Friedrich et al, Phys. Rev. Lett. 109 138102, 2012), a minimal model swimmer 
was proposed that propels itself at low Reynolds numbers by a revolving motion of a pair of spheres. 
The motion of the two spheres can synchronize by virtue of a hydrodynamic coupling that depends on 
the motion of the swimmer, but is rather independent of direct hydrodynamic interactions. This novel 
synchronization mechanism could account for the synchronization of a pair of flagella, e.g. in the green 
algae Chlamydomonas. Here, we discuss in detail how swimming and synchronization depend on the 
geometry of the model swimmer and compute the swimmer design for optimal synchronization. Our 
analysis highlights the role of broken symmetries for swimming and synchronization. 
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I. INTRODUCTION 



Hydrodynamics at the scale of a single cell can be counter-intuitive, as inertia is negligible and 
propulsion by thrust becomes impossible [1-3]. In this regime of low Reynolds numbers, the 
hydrodynamic equations are symmetric under time-reversal and net propulsion requires a swim- 
ming stroke that breaks time-reversal symmetry [4] . In pioneering work, Taylor demonstrated 
theoretically that self -propulsion by purely viscous forces is indeed possible when he studied a 
minimal system of waving elastic sheets that abstracts from the undulatory bending waves of 
the flagella of sperm cells [5]. Since then, simple model swimmers became an indispensable 
tool to explore basic principles of swimming at low Reynolds numbers, both theoretically [6-9] 
and experimentally [10-12]. 

Already in his early work [5], Taylor addressed also the striking biological phenomenon of syn- 
chronization among beating flagella. Eukaryotic flagella are long, slender appendages of many 
cells [13, 14] that beat actively due to the collective activity of molecular motors inside [15]. 
The active periodic bending waves of several flagella can synchronize, as observed e.g. for the 
two flagella of a pair of nearby sperm cells, [16, 17], or among the hundreds of short flagella 
on ciliated epithelia in mammalian airways [18]. Recently, the green algae Chlamydomonas, 
a unicellular organisms that swims with a pair of flagella like a breast- swimmer [19], emerged 
as an experimental model system for flagellar synchronization [20-23]. In Chlamydomonas, 
synchronized flagellar beating is required for swimming along a straight path [19, 21]. Coor- 
dination of flagellar beating is physiologically important for cellular propulsion also in other 
systems [13, 24], as well as for efficient fluid transport [25, 26]. Taylor suspected that flagel- 
lar synchronization is a hydrodynamic effect [5], although a definite proof is pending to date 
[27]. Again, simple model systems served as a proof of principle that synchronization by hy- 
drodynamic forces is possible, both in theory [28-33] and experiment [34-38]. These model 
systems demonstrated equally clearly the need for non-reversible phase dynamics and thus bro- 
ken symmetries [39]. In these models, symmetries were broken either by wall effects [28], by 
additional elastic degrees of freedom [29, 30], or by evoking phase-dependent driving forces 
[33, 36]. Here, we study theoretically a novel class of simple swimmers that are inspired by the 
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propulsion of Chlamydomonas cells and comprise a pair of driven spheres instead of two flag- 
ella (a special case was reported in [40], see also [41]). The two driven spheres of the swimmer 
can synchronize their motion by hydrodynamic forces in a way that is independent of hydro- 
dynamic interactions between them, but crucially depends on the swimmer's ability to move: 
If the two driven spheres are initially out of phase, an imbalance of local torques causes the 
whole swimmer to rotate (and to translate). This motion imparts different hydrodynamic fric- 
tion forces on the two spheres, causing one of them to slow down its orbit, and the other one to 
speed up, which ultimately restores synchrony. Like swimming, synchronization requires that 
time-reversal symmetry is broken, which puts certain constraints on the design of the swimmer. 
Our model swimmer exemplifies a generic mechanism for hydrodynamic synchronization, 
which differs in a fundamental way from other mechanisms described previously [28-34, 37] 
that depended on hydrodynamic interactions. This hydrodynamic synchronization mechanism 
suggests a specific physical mechanism for flagellar synchronization in Chlamydomonas. Fur- 
ther, we anticipate that this synchronization mechanism could be exploited for the design of 
novel man-made micro- swimmers. 

Outline. In section 2, we specify the geometry of the model swimmer as well as its equation of 
motion, which couples the phase dynamics of its swimming strokes and the swimming motion. 
Following [40], this equation of motion is derived as a force balance between internal active 
driving forces and generalized hydrodynamic friction forces using the framework of Lagrangian 
mechanics for dissipative systems. In section 3, we discuss self-propulsion and self-rotation of 
the swimmer, while section 4 addresses the phenomenon of phase synchronization of its two 
driven spheres in relation to the swimmer's ability to move. 

H. A MINIMAL THREE-SPHERE SWIMMER FOR FLAGELLAR SYNCHRONIZATION 

Inspired by a biological micro- swimmer, the biflagellate green algae Chlamydomonas, we pro- 
posed a minimal swimmer that retains the basic symmetries of a Chlamydomonas cell and is 
likewise propelled by two 'flagellar actuators' [40]. This model swimmer is shown in figure 
1 and consists of three spheres, which are connected by a frictionless planar scaffold. Two of 
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the three spheres can move along circular trajectories being connected by lever arms of length 
R to the corners s*, % = 1, 2 of an isosceles triangle. These two spheres have equal radius a 
and mimic the two flagella of Chlamydomonas. The third sphere of radius b plays a role anal- 
ogous to the body of this cell. For simplicity, we assume that the swimmer is planar and that 
the positions of the three spheres all lie in the rry-plane. The position and spatial orienta- 
tion of each sphere is thus characterized by just two coordinates r; = (x i: y i: 0) with respect 
to the rryz-laboratory frame as well as a single rotation angle a, for rotations around the z- 
axis. For brevity, x = x 3 , y = y 3 , a = a 3 . We explicitly note a material frame for the third 
sphere, which consists of the orthonormal vectors ei = (cos a, sin a, 0), e 3 = (0,0,1) and 
e 2 = e 3 x ei. The geometry of the swimmer is characterized by the dimensions / and h of 
the triangular scaffold such that s, = r 3 + (— l) l Zei + he 2 , i — 1, 2, see again figure 1. The 
circular orbits of the two flagellar spheres are parametrized by respective flagellar phase angles 
(fi with r; = Sj + R(— snivel + cosy9;e 2 ). Here, the flagellar phase angles are defined as 
continuous variables such that the number of full cycles of the ith sphere after a time t is given 
by \_(ifi(t) — y?j(0))/(27r)J, % — 1,2. Likewise, the number of full rotations of the ith sphere 
with respect to the laboratory frame is given by \ ioLiit) — aj(0))/ (2%)\ ,i — 1, ... ,3, where we 
have the identity aij = a + ipi, i — 1, 2, for the first and second sphere. Throughout the text, 
positive angular velocities <fii indicate counter-clockwise rotations (when viewed along — e 3 ), 
while clockwise rotations correspond to negative angular velocities of the respective sphere. 



A. Low Reynolds number hydrodynamics 

For a swimming microorganism, inertial forces are usually negligible in comparison to viscous 
forces, which is reflected by a small Reynolds number [2, 3] 

Re=^. (0.1) 
V 

Here, I is a length-scale of the swimmer, v its speed, and 77 and p the viscosity and density of 
the surrounding fluid. As an example, we estimate Re ~ 10~ 3 for the unicellular green algae 
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Chlamydomonas using v ~ 100 /im/s and / ~ 10/im [19]. 

In the following, we will always assume that inertial forces are also negligible for our model 
swimmer. The hydrodynamics at low Reynolds numbers is then governed by the Stokes equa- 
tion, Vp + r^V 2 u = 0, where p, u, 77 denote the pressure field, flow field and viscosity of the 
fluid, respectively [2]. Importantly, this equation is linear, which implies a linear relationship 
between velocities and resultant forces. For example, a single sphere of radius a that moves 
with translational velocity v through a viscous fluid will exert a hydrodynamic friction force 
F = 7V on the fluid, where 7 = Giir/a is the Stokes friction coefficient for translational mo- 
tion [2]. Likewise, a rotation with angular velocity u> induces a hydrodynamic friction torque 
T' = 7 r a; with rotational friction coefficient 7 r = 8rrr]a 3 . In the following, primed torques will 
always denotes torques relative to the center r of the respective sphere. 

For an ensemble of N spheres with positions rj and radii a i5 the linearity of the Stokes equation 
implies a similar linear relationship [2, 3, 42] 

Po = r • q , (112) 

where q is a 6iV-vector that combines the translational and rotational velocity com- 
ponents of the N spheres, qo = (vi x ,vi y ,vi z ,0Ji x ,0Jiy,0Ji z , . . . ,un z ), while the QN- 
vector P combines the components of the resultant hydrodynamic friction forces, P = 
(F lx , F ly , F lz , T[ x , T[ y , T[ z , T' Nz ). Finally, T denotes the 6N x 6N grand friction ma- 
trix. In the limit of large separation = | r» — r^- 1 3> a { between the spheres, the friction matrix 
becomes diagonal; the diagonal is then populated by the common Stokes friction coefficients 
for translational and rotational motion of single spheres. Generally, off-diagonal entries of r 
characterize hydrodynamic interactions between the spheres: Any movement of the ith sphere 
will set the surrounding fluid in motion. If the jth is held fixed in the vicinity, it will slow down 
this fluid flow and exert a corresponding friction force on the fluid. 

The friction matrix T is always symmetric as a consequence of the Lorentz reciprocal theorem 
[2]. Furthermore, T is positive definite, as the hydrodynamic dissipation rate lZ h = ■ P 
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can never become negative [2] 

K h = cil • r • q > 0. (II.3) 

The friction matrix T can be computed to arbitrary precision in aj/r^ by using e.g. the method 
of reflections [42, 43]. The leading order contribution is provided by pairwise interactions 
between the spheres and amounts to the Rodner-Prager approximation, which generalizes the 
popular Oseen tensor, see appendix A. 

B. Over-damped Lagrangian mechanics 

We employ the framework of Lagrangian mechanics for dissipative systems [44, 45] to derive 
equations of motions of our model swimmer from simple force balances [40]. The swimmer is 
characterized by just 5 degrees of freedom if its internal constraints are taken into account: The 
two phase angles <f 1 and y? 2 describe the orbiting of the two flagellar spheres, while x, y and a 
characterize translations and rotations of the swimmer in the plane of swimming. We introduce 
a 5-component vector of generalized coordinates 

q = (x,y,a,(p 1 ,(f 2 ) T . (II.4) 

The full position vector q for the N = 3 spheres can now be uniquely expressed as a function 
of the generalized coordinates as q = qo(q)- For each degree of freedom, we define a general- 
ized hydrodynamic friction force Ph,i conjugate to this degree of freedom such that the rate of 
hydrodynamic dissipation can be equivalently written as 

K h = P h ■ q = P h ,!X + P h , 2 y + P h ,sa + P h ,m + P h ,5^2, (H.5) 

where the 5-vector P h = (Ph,i, ■ ■ ■ ,Ph,b) is simply given by P h = L T • P with a 9 x 5 
transformation matrix Lij = dq ,i/dqj. The generalized hydrodynamic friction forces P hyi 
relate linearly to the generalized velocities <jj as P h)i = Yfj=i ^h,ijQj- Here, we have introduced 
the generalized 5x5 hydrodynamic friction matrix T h = L T • T ■ L, which is just a contracted 
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version of the grand friction matrix r . Note that q = L ■ q. 

We further include internal dissipation related to the motion of the two driven spheres with 
associated dissipation rate 1Z K = /t(y5i 2 + y?2 2 )- The total dissipation rate is now given by 
1Z = 1Z h + 1Z K and defines a Rayleigh dissipation function for our mechanical system. We may 
write 

11 = q T • r • q, (II.6) 

where for notational convenience T is defined as T = T h +T K with a 5 x 5 internal friction matrix 
T K , whose only non-zero entries are T K j44 = T K 55 = k. The Rayleigh dissipation function is 
related to the generalized friction forces P by 2Pj = dlZ/dqi, i = 1, . . . , 5. Explicitly, these 
generalized friction forces read 

Pi = Pi,z + F 2)X + F 3tX , P 2 = F lt y + F 2tV + F 3>y , (H.7) 
P 3 = [(n - r 3 ) xF 1 + (r 2 - r 3 ) x F 2 ] z + T[ z + T' 2z + T^, (II.8) 
P 3+z = ¥ r {dv 3 /d^)+T' hZ + K^ j = 1,2. (11.9) 

Note that Pi and P 2 simply denote the x- and ^/-component of the total hydrodynamic friction 
force exerted by the swimmer on the fluid. Likewise, P 3 denotes the ^-component of the total 
torque exerted by the swimmer on the fluid (with respect to r 3 ). For free swimming, the swim- 
mer is free from external forces and torques, which implies a force balance, P 1 = P 2 = as 
well as a torque balance P 3 = 0. 

The two flagellar spheres are actuated by generalized active driving forces Q A and Q 5 (with 
have units of a torque). For simplicity, we assume these generalized forces to be constant, 
independent of the flagellar phase, Q A = mi and Q 5 = m 2 . The active driving forces can be 
re-written as potential forces: The energy for active swimming is provided by a fuel reservoir 
with internal energy U such that — U = 1Z. The simplest choice U = —m\ipi — m 2 tp 2 defines 
generalized potential forces Qi = —dU/dqi, i = 1, ... ,5, and recovers the above definition 
with phase-independent Q A and Q 5 . A case of phase-dependent driving forces was recently 
studied numerically in [41]. If inertial forces can be neglected as in our case and the swimmer 
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is free from external forces and torques, we obtain a balance of generalized forces, Pi = Qi, 
i = 1, . . . , 5, where Qi = Q 2 = and Q% = 0. In matrix notation, this force balance reads for 
this case for free swimming 

(0,0,0,m 1 ,m 2 ) T = r-q. (11.10) 

C. Equations of motion 

Free swimming. The balance of generalized forces stated in eq. (11.10) enables us to self- 
consistently solve for the dynamics q of the swimmer. It is convenient to introduce a block 
matrix notation for the generalized friction matrix 

(k c\ 

T= ; (11.11) 

\c T nj 

here K, C and Q have dimension 3 x 3, 3 x 2 and 2x2, respectively. Similarly, we write 
q = (X, <&) with X = (x, y, a) and $ = (</?i, (p 2 ). From equation (11.10), we find 

X = -K 1 C$, (11.12) 
$ = (fi-C T KC)- 1 (m 1 ,m 2 ) T . (11.13) 

The first equation (11.12) states that the motion X of the swimmer is fully determined by its 
phase dynamics The phase dynamics $ in turn depends only on but not on the position 
and orientation of the swimmer as encoded by X. In fact, the friction matrix T is independent 
of the position of the swimmer. Furthermore, a rotation of the swimmer around e 3 with rotation 
matrix D transforms T according to K — > D _1 KD, C — > D _1 C, ft — > ft. In particular, 
ft — C T KC does not change under such a rotation. 

Clamped third sphere. If the third sphere is constrained from translating and rotating, we natu- 
rally have X = 0. In this case, a generalized constraining force C$ must be applied at r 3 . The 
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phase dynamics is then given by 



(m 1 ,m 2 ) T . 



(11.14) 



Both equation (11.13) and equation (11.14) describe a pair of coupled phase oscillators [46] and 
can be recast in the form 



Importantly, the coupling functions Hi are different in the case of free- swimming and the case of 
a clamped third sphere. In the case, where the third sphere is clamped, the coupling functions 
Hi comprise only direct hydrodynamic interactions between the spheres. In the case of free 
swimming, the coupling functions Hi implicitly account also for the motion of the swimmer 
and the associated hydrodynamic friction forces acting on the two flagellar spheres. 

III. SWIMMING BY HYDRODYNAMIC INTERACTIONS 
A. Swimming forward or backward 

We study the case of exactly opposite driving torques m 2 = —mi, which gives rise to a counter- 
rotation motion of the first and second sphere, similar to the mirror-symmetric beat of the two 
flagella of Chlamydomonas. In this case, there exists a limit cycle of the phase dynamics with 
5 = <f i + (f 2 = 0, which corresponds to in-phase motion of the two spheres. In section IV, we 
show that this limit cycle is stable for ui h < and unstable for cu h > 0. Here, we introduced 
the (signed) angular frequency cu = mi/ k. The corresponding period is T = 27r/w . 
To gain analytical insight, we resort to a limit of small spheres and short lever arms lengths. We 
introduce a small expansion parameter e = a/l and assume b/l and R/l to be of order unity. 
Details of the calculation can be found in appendix B. 

For mirror-symmetric propulsion with 5 = 0, the swimmer will move along a straight path with 
velocity v = r 3 • e 2 and zero rotation a — 0. 



Kifii = mi + Hi{ipx,ip 2 ), 



i = 



1,2. 



(11.15) 
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Remarkably, even for zero lever arm length R = 0, when the first and second sphere spin only 
around their respective centers, the swimmer is able to move forward by harnessing hydrody- 
namic interactions that couple to rotational motion, see also [47] for a similar swimmer design. 
For the swimming speed, we find in this case 

i , a ( a 2 2abl \ . 4 . 

fcr = ^ where v ° = u ° a 2a-Tb \w + w + ) + ° {£ ] - (ml) 

Note that net propulsion is possible even in the absence of the third sphere, 6 = 0, which 
provides a simple design of a two-sphere-swimmer. 

In the general case of a finite lever arm length R > 0, the swimmer will periodically wiggle 
back-&-forth to leading order in e 

v = u R^- sin^i + 0(e 2 ). (111.2) 
2a + b * v ' 

Net propulsion is a higher-order effect with time-averaged velocity (v) = 0(e 3 ) that requires 
hydrodynamic interactions [2, 45]. In the limit of small lever arm lengths, R<^l,we find 

(v) -v ~R 2 + 0(R 3 ), (III.3) 

where v = lim R _ >0 (f). Eq. (III.3) exemplifies the basic fact that the net swimming speed, 
which a swimmer can achieve by periodic body shape changes, generically scales with the 
square of the amplitude of its swimming stroke. This "quadratic law of low Reynolds number 
propulsion" holds for other model swimmers [6, 7, 9] and was formally proven by Shapere and 
Wilczek [4] by generalizing Purcell's scallop theorem [1]. 

While non-zero net propulsion is a generic consequence of a non-reciprocal swimming stroke 
that breaks time-reversal symmetry, the actual sign of the velocity can depend on the geometry 
of the swimmer [48]. In our case, (v) can switch sign as a function of b/a and h/l, see figure 2. 
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As a consequence, our swimmer may either swim forward or backward for w > 0, 

<„> = ,o - ab "° R2 ^ %f+ a m 2 b + for h = h (HI.4) 

(v) =v + abu R 2 J°^-^ 2 + 0(e 4 ) for h = 0. (HI.5) 

8(2a + byV 

This surprising behavior can be related to two opposing mechanisms that contribute to net 
propulsion in the case of in-phase beating with 5 = 0: The instantaneous swimming speed 
v — f/g can be computed from the friction coefficient g associated with towing the swimmer 
along the e 2 -direction, and a force — fe 2 that equals the force that had to be applied at r 3 to 
prevent the swimmer from moving. (If a — 0, g — K 2 2 and (0, /, 0) T = — C<&.) The force 
/ oscillates during each swimming stroke, / = 12ttt] a Ru sin (fi + 0(e 3 ), which is tightly 
linked to the back-&-forth motion of the swimmer as characterized by equation (III. 2). Now, 
hydrodynamic interactions between the spheres result in higher-order corrections to / and im- 
ply in particular (/) > for cu > 0. This can be seen as follows: When the spheres are close 
together (a configuration for which / < if lu > 0), hydrodynamic interactions are strongest 
and reduce the magnitude of /. Averaging over one phase cycle thus leaves a positive net con- 
tribution. This is the first mechanism, which contributes to forward swimming for uj > 0. At 
the same time, the friction coefficient g depends on the phase (pi. it is largest, when the spheres 
are furthest apart, and smallest, when they are close together. The oscillations in the friction 
coefficient g can resonate with the oscillations of the force / and contribute to backward swim- 
ming. Depending on the geometry of the swimmer, either mechanism may dominate resulting 
in either (v) > or (v) < 0, see figure 2. 

B. Asynchronous beating results in rotational motion 

In the case of a phase-difference between the two driven spheres, 5 = <fi + v?2 0, mirror- 
symmetry with respect to the central axis of the swimmer is broken and the swimmer can rotate 
with a 0, where a — a s for brevity. These rotations can feed back on the phase dynam- 
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ics of the two spheres and cause them to synchronize as discussed in the next section. In this 
subsection, we want to elucidate how asynchrony causes the swimmer to rotate, without further 
complications arising from the feedback on the phase dynamics. To this aim, we will temporar- 
ily prescribe the phase dynamics as ip 1 = u t, (f 2 = —u>ot + 5. We then find that the swimmer 
"rocks" in a periodic fashion, 



Figure 3 further shows the net rotation rate, averaged over one beat cycle. In the special case, 
where I — h, b — a, we have 



A qualitatively similar behavior is found for sideways translation motion with velocity r 3 • ei. 
Role of symmetries. The net rotation rate (a) is always an odd function of the phase difference 
5. This behavior is a direct consequence of the symmetries of the swimmer: A reflection at 
the central symmetry axis of the swimmer maps the swimmer onto a copy of itself, but with 
phase difference 5 — > —5. Since a reflection changes the sign of a, the assertion follows. As a 
corollary, the swimmer will show zero net rotation for 5 = it. In a similar way, (a) is an odd 
function also of the position h of the third sphere. To see this, consider a rotation that rotates 
the swimmer in the rry-plane by an angle n around the point (si + s 2 )/2. This rotation maps the 
swimmer onto a swimmer for which both h and u have the opposite sign. Of course, a does 
not change under this operation, which implies that (a) is an even function of u h. Since uj 
sets the time-scale of the problem, we have a ~ u and the assertion follows. As a particular 
case, the swimmer will show zero net rotation if h — 0. In the next section, we will encounter 
a similar behavior for the synchronization strength. 



a = 



u Rsin(S/2) 



[(2a + 6)/ cos(v?i - 5/2) - bh sin(v?i - 5/2)} + 0(e 2 ). (III.6) 



(2a + b)P + bh 2 
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IV. SYNCHRONIZATION DEPENDS ON MOTION 

In our model, the phase velocity of each of the two driven spheres, (pi, is determined by a 
balance of an active driving force and a generalized hydrodynamic friction force -Pj+3, 
% = 1,2. The hydrodynamic friction forces P; are affected by the instantaneous motion of the 
swimmer, which indirectly results in a coupling between the two driven spheres, see section 
IIC. To characterize the resultant phase dynamics, we resort to a stroboscopic description at 
times t n that mark the completion of a full cycle of the first sphere, such that ipi(t n ) = 2nn 
with integer n. We introduce a Poincare return map A(5) that measures the change of the phase 
difference 5 = <f± + y?2 after a full cycle of the first sphere, A(5(t )) = 5(ti) — 5(t ). To leading 
order, we find [40] 



Here, A = — dA/d5\s=o denotes a dimensionless synchronization strength [49]: The limit cycle 
of in-phase synchronized phase dynamics with 5 = will be exponentially stable with Lya- 
punov exponent ln(l — A) whenever A > 0. For free swimming, we find A > provided 

uoh < 



see also figure 4D. Expression (IV.2) for the synchronization strength A has a single global 
maximum as a function of the size b and position h of the third sphere, which is characterized 
by b = a and h — l. Thus, the swimmer design shown in figure 1 is optimal for synchronization 
in the case of free swimming. Interestingly, hydrodynamic interactions contribute only to higher 
order to the synchronization strength A [40] . Instead, synchronization strongly depends on the 
swimmers ability to move. For our particular swimmer design, strong synchronization with A = 
0(s 5 ) requires both translational and rotational motion. As a means of illustration, consider the 
case, where translation of the swimmer is constrained by an external force, while the swimmer 



A(5) 



(IV.l) 




(IV.2) 
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can still rotate freely. We then find A = 0(e & ), more specifically 

A = -signM*) ^ [l7 + 5(V0 2 -8(a/i?) 2 (3 + (V0 2 )] ( J^f^ ) ' + 0(e r ). 

(IV.3) 

In this case, hydrodynamic interactions contribute to leading order. Similarly, if the swimmer 
can only translate, but not rotate, we obtain 

, lN 30afo|% l(l 2 -h 2 ) ( 3naR 2 \ 2 ml 7n 

Finally, if the swimmer is fully clamped and can neither translate nor rotate, we have A = 0(s 8 ) . 
Role of symmetries. For any dynamics of the swimmer, we can realize the time -reversed dy- 
namics by simply reversing the sign of the active driving forces, m x — > —mi and m 2 — > —m 2 . 
Under this mapping, a stable limit cycle of synchronized phase dynamics will become unstable, 
and vice versa. Generally, synchronization requires that the phase dynamics is non-reversible 
[39], which implies that the time -reversed dynamics should not be equivalent to the original 
dynamics. For our model swimmer with mi = — m 2 , the phase dynamics becomes reversible 
if h = 0. As a proof, consider again the rotation that rotates the swimmer in the rry-plane by 
an angle n around the point (si + s 2 )/2. The rotated swimmer is equivalent to a swimmer with 
mi — > m 2 = —mi, m 2 — > mi — —m 2 and h — > —h. On the other hand, a mere rotation 
of the swimmer does not change its synchronization behavior. We conclude that the phase dy- 
namics of the original swimmer is equivalent to the time -reversed dynamics of a swimmer with 
h — > —h. For h — 0, the phase dynamics is thus reversible and all orbits of the phase dynamics 
are neutrally stable. In particular, A = for h = 0. As a corollary, the synchronization strength 
A also vanishes in the absence of the third sphere, since A is independent of h in this case. The 
situation is different for a swimmer that cannot translate, but still rotate around r 3 . The very 
fact that the swimmer is acted upon by constraining forces to ensure r 3 = breaks reversibility 
and results in stable synchronization even in the absence of the third sphere, provided cu h < 0, 
see eq. (IV.3). 



14 



V. DISCUSSION 



In this manuscript, we presented a simple model swimmer to demonstrate a generic mechanism 
for hydrodynamic synchronization. This synchronization mechanism depends on the bidirec- 
tional coupling between the phase dynamics of two oscillators that propel the swimmer and 
the resultant swimming motion. Including the positional and orientational degrees of free- 
dom of the swimmer is sufficient to break the reversibility of the phase dynamics, which is a 
prerequisite for synchronization [39]. Restraining the swimmer's ability to move reduces the 
synchronization strength in a characteristic manner, resulting in different scaling behaviors of 
the synchronization strength for different constraining conditions. From our analysis, we can 
furthermore deduce the swimmer design that is optimal for synchronization. The existence of 
an optimal design can be intuitively understood as follows: In the absence of the third sphere 
(which represents a passive drag center), the swimmer attains an additional symmetry, which 
rules out synchronization. On the other hand, a third sphere that is very large in comparison to 
the two actively driven spheres would attenuate the swimming motion of the swimmer and thus 
reduce its synchronization strength. In fact, we show in a limit of small spheres that synchro- 
nization works best if all three spheres are of same size (b = a). A similar reasoning applies 
also to the position of the third sphere. Synchronization is ruled out by symmetry in a swimmer 
design for which the third sphere is collinear with the rotation centers of the two driven spheres. 
For a large distance between the third sphere and these rotation centers, however, rotations of the 
swimmer are significantly reduced, which in turn implies a reduced synchronization strength. 
Correspondingly, we find that the synchronization strength is maximal for an intermediate value 
of this distance (h = I). 

In addition to synchronization, our model swimmer allows us to illustrate generic features of 
self-propulsion at low Reynolds numbers. For example, our swimmer displays characteristic 
back-&-forth wiggling on top of net propulsion, which represents a hallmark of self -propulsion 
at low Reynolds numbers by cyclic swimming strokes. Unlike net propulsion, this periodic 
motion requires neither non-reversible swimming strokes nor hydrodynamic interactions. This 
periodic motion is crucial for the synchronization mechanism studied here. A wiggling back- 
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&-forth motion is not only predicted for model swimmers, but actually observed for biological 
micro-swimmers that use cyclic swimming strokes. For example, sperm cells display a char- 
acteristic wiggling motion of the sperm head, being propelled by undulatory bending waves of 
their eukaryotic flagellum [50, 51]. Another class of biological swimmers, however, such as 
Paramecium generate slip flows at their boundary to swim at a steady velocity without apparent 
body changes [13]. Our simple three-sphere swimmer incorporates both propulsion primitives 
as limit cases. 

In conclusion, we have shown that swimming and synchronization can strongly depend on 
each other. Both net propulsion and net synchronization are a generic consequence of broken 
symmetries. The actual sign of the swimming speed or synchronization strength, however, will 
depend on the geometry of the swimmer in a non-generic way. 

Appendix A: Rodner-Prager approximation of hydrodynamic interactions 

In this appendix, we show how the grand hydrodynamic friction matrix r , can be computed 
to leading order approximation for an ensemble of iV spheres of different radii, see [42] for 
the case of equally-sized spheres. The inverse of r , the grand mobility matrix M = Tq 1 , 
describes the linear dependence between the translational and rotational velocities Vj and Wj of 
the spheres, % — 1, . . . , N, and the forces Fj and torques T; exerted by the spheres on the fluid, 

N 

v, = XX-F,- + mJT,-, (A.l) 

N 

Here, we have decomposed M into 3x3 blocks as M = (/it", fJ-ij)i,j=i,...,N to show 

the different types of couplings explicitly. With the compact notation from section II A, we 
could equivalently write q = M • P . 

The method of reflections can be used to construct approximations of arbitrary precision of M 
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[42, 43]. The leading order approximation for fj,^ is known as the Rodner-Prager approximation 
and treats a sphere with index j as a test particle that is convected by the fluid flow induced by 
another sphere with index % that moves under the influence of a force (or torque). Higher-order 
corrections would for example describe how the presence of the jth sphere changes the mobility 
of the ith sphere. 

We denote the sphere's positions by r; and their radii by i = 1, . . . , N. The flow field u* 
induced by a single sphere of radius a« translating under the influence of an external force Fj 
with a constant velocity = Fi/%, ^ = 6-irrjai, in a viscous fluid of viscosity rj reads 

u\(r,t) = A(r- ri ) Vl , A(r) = ^(l + r®r) + ^^) 3 (l-3rg)r). (A.3) 

Here r = |r — | , r = (r — r^/r. If the sphere is subject to an external torque Tj and rotates 
with rotational velocity u>j = Tj/7[, Yi — Sirr/af, then the induced flow field is 

<(r)= (^) 3 «xr. (A.4) 

Faxen's theorem provides the translational and rotational velocities of a second sphere with 
label j that is free from external forces and torques and becomes convected by the flow Uj 

Vj = £*u(r)| r=r „ C* = l+ l -a)V\ (A.5) 
^.=/7u(r)| r=rj , r = ^Vx. (A.6) 

Combining the explicit expressions for the flow fields induced by the ith sphere given by 
eqns. (A.3),(A.4) and Faxen's laws provides leading order approximations for the cross- 
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mobilities jj, h [42, 43] 



H% = ^r (l + *M + - + 0(e*) (A.7) 

1 



^ = -16^( 1 -^-) +^ 6 ) (A. 8) 

Mg = ^ = ^^x +(9(, 5 ). (A.9) 



Here, we introduced a small expansion parameter e and assumed that all sphere radii are small 
of order e compared to the distances r^. for the rotational-translational coupling. The self- 
mobilities are perturbed only to higher order [42, 43] 

H* = ±l + 0(e% i%=± : l + 0(e% K' = Oie 7 ), (A. 10) 

Appendix B: A limit of small spheres and lever arms 

To gain analytical insight into the equation of motion, we resort to a limit of a small beat am- 
plitude characterized by small lever arms and small spheres. We introduce a small expansion 
parameter e = a/l and assume both b/l and R/l to be of order unity. We expand the force bal- 
ance given by eqn. (11.10) as a power series in e making use of the block matrix decomposition 
(11.11) 

/o\ 




mi 
\m 2 J 



£K« +E 2 K^ +... £ 2 C (2) +£ 3 C (3) + \ / £ x(l) +£ 2 X (2) + 
£ 2 C (2),T +£ 3 C (3),T + _ _ _ ^(0) +e 3^(3) + _ _ I I $(0) +£ 3$(3) + _ _ _ 



Note that and are invertible. We can now iteratively solve for the series coefficients 
of X and Hydrodynamic interactions between the three spheres contribute to leading order 



18 



only to K( 2 ), C( 3 ) andfi^. 
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FIG. 1: A simple three-sphere swimmer. A. The swimmer consists of three spheres connected by a 
frictionless, planar scaffold. The first and second sphere (red, located at ri and r-2) can move along 
circular orbits, thus propelling the swimmer. The swimmer is characterized by 5 degrees of freedom: 
two phase angles tpi and ip% that parametrize the circular orbits of the first and second sphere, as well 
as the planar position (x, y) and orientation angle a of the third sphere (blue, located at r3). B. Flow 
fields induced by the swimming motion of the swimmer when immersed in a viscous fluid. Parameters: 
a = b = 0.1 1, h = 2R = I, mi = n = r]l 3 . 
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FIG. 2: Forward or backward swimming for different swimmer geometries. For synchronized phase 
dynamics with S = ip% + tp2 = 0, the swimmer will move along a straight path with velocity v = ■ e2. 
The net swimming velocity, averaged over a beat cycle, depends on the radius b of the third sphere (A) 
as well as on its position h (B). Solid lines indicate the analytical approximation, while the dashed line 
was determined from numerical integration of the equation of motion. Parameters: As in figure 1, unless 
indicated otherwise. 
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FIG. 3: Asynchronous beating causes rotation. For this figure only, we us a prescribed phase dynamics 
with ipi = ujQt and if2 = —coot + 5 to separate the swimming motion from its feedback on the phase 
dynamics. A. Swimming trajectory T3(t) for constant phase difference 5 = ir/2. The position r^itn) 
and orientation e2(t n ) after each full cycle of the first sphere, ipi(t n ) = 2nn, is indicated as dot and 
arrow, respectively. B. The orientation angle a = a^(t) oscillates around a linear trend, indicating non- 
zero net rotation after each full cycle. C. Different swimmer designs studied in panel D: top (red): the 
design from figure 1, middle (blue): h = 0, bottom (green): without the third sphere, b = 0. D. For 
the swimmer design from figure 1, the net rotation rate (a) depends on the phase difference 5 as well as 
on the radius b and position h of the third sphere (red solid line: analytical approximation, red dashed: 
numerical result). Additionally, we show numerical results for the case of scaled-down spheres and lever 
arms with a = b = el, R = 5el using e = 0.01 (red dotted curve). The net rotation rate vanishes due to 
symmetry reasons for b = (blue curve) and h = (green curve). Parameters: As in figure 1-2, unless 
indicated otherwise; panel D, middle and right: 5 = tr/2. 
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FIG. 4: Phase-synchronization of a free-moving swimmer. A. Swimming trajectory rs(t) for a swimmer 
with loq < and initial phase difference 5(0) = ir/2. The position r^(t n ) and orientation e2(t n ) after 
each full cycle of the first sphere, <pi(t n ) = —2-kti, is indicated as dot and arrow, respectively. B. The 
orientation angle a = 03(f) oscillates with an amplitude that decays in time as the two driven spheres 
synchronize in-phase. C. The phase-difference S = <~p\ + decays towards zero, indicating in-phase 
synchronization. D. Stable in-phase synchronization is characterized by a synchronization strength A > 
(solid line: analytical approximation eq. (IV.2), dashed line: numerical result). Here, A is defined as 
A = — dA(S) / dS\g =0 , where A(<5(irj)) = S(t±) — 5 (to) measures the change in phase difference after one 
full cycle of the first sphere. In a limit of small spheres, e<1, synchronization strength A is maximal 
for a swimmer design with b = a, h = I as indicated. In the absence of the third sphere, 6 = 0, the 
synchronization strength vanishes, since the dynamics of the swimmer becomes reversible in this case. 
Similarly, one can derive the symmetry relation ln[l + A(— h)] = — ln[l + X(h)], which implies that the 
in-phase synchronized state with 6 = changes its stability if h reverses sign. Parameters: As in figure 
1-3, except that the sense of rotation of the first and second sphere is reversed, —mi = ni2 = K = rjl 3 . 
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